Self-organization of traffic jams in cities: effects of stochastic dynamics and signal 
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We propose a new cellular automata model for vehicular traffic in cities by combining (and appro- 
priately modifying) ideas borrowed from the Biham-Middleton-Levine (BML) model of city traffic 
and the Nagel-Schreckenberg (NS) model of highway traffic. We demonstrate a phase transition 
from the "free-flowing" dynamical phase to the completely "jammed" phase at a vehicle density 
which depends on the time periods of the synchronized signals and the separation between them. 
The intrinsic stochasticity of the dynamics, which triggers the onset of jamming, is similar to that 
in the NS model, while the phenomenon of complete jamming through self-organization as well as 
the final jammed configurations are similar to those in the BML model. Using our new model, we 
have made the first investigation of the time-dependence of the average speeds of the cars in the 
"free- flowing" phase as well as the dependence of flux and jamming on the time period of the signals. 
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Over the last half century various concepts and tech- 
niques of fluid dynamics and statistical mechanics have 
been successfully applied to understand several funda- 
mental aspects of vehicular traffic flow ||^j)> The 
"particle- hopping" models [|][5| of vehicular traffic are 
usually formulated using the language of cellular au- 
tomata (CA) Q. These models are closely related to 
some of the microscopic models of driven systems of in- 
teracting particles, which are of current interest in non- 
equilibrium statistical mechanics 0. 

A one-dimensional CA model of highway traffic and a 
two-dimensional CA model of city traffic were developed 
independently by Nagel and Schreckenberg (NS) M and 
Biham, Middleton and Levine (BML) ||, respectively. 
Highway traffic becomes gradually more and more con- 
gested in the NS model with the increase of density. Traf- 
fic jams appear in the NS model because of the intrinsic 
stochasticity of the dynamics H but no jam persists for 
ever. On the other hand, a first order phase transition 
takes place in the BML model at a finite non-vanishing 
density, where the average velocity of the vehicles van- 
ishes discontinuously signalling complete jamming. In 
the BML model, the randomness, which is crucial for the 
jamming, arises only from the random initial conditions, 
as the dynamical rule for the movement of the vehicles is 
fully deterministic ||. 

If each unit of discrete time interval in the BML model 
is interpreted as the time for which the traffic lights re- 
main green (or red) before switching red (or green) si- 
multaneously in a synchronized manner, then, over that 
time scale each vehicle, which faces a green signal, gets an 
opportunity to move from one crossing to the next. The 
generalization of the BML model that we propose here 
is, to our knowledge, the first attempt to describe explic- 
itly the forward movement of the vehicles over smaller 
distances during shorter time intervals. We achieve this 
generalization by following the prescriptions of the NS 
model not only for describing the positions, speeds, ac- 
celerations and decelerations of the vehicles |j but also 
for taking into account the interactions among the ve- 
hicles moving along the same lane of a street. We also 



modify some of the prescriptions of the BML model ap- 
propriately to take into account the signal-vehicle inter- 
actions and the interactions between vehicles approach- 
ing a crossing along different streets. 

Our main aim is to demonstrate that a phase transi- 
tion from the " free-flowing" dynamical phase to the com- 
pletely "jammed" phase takes place in our generalized 
model; the intrinsic stochasticity of the dynamics, which 
triggers the onset of jamming, is similar to that in the 
NS model, while the phenomenon of complete jamming 
through self-organization as well as the final jammed con- 
figurations are similar to those in the BML model. 

In the BML model a square lattice represents the net- 
work of the streets. All the streets parallel to the X- 
direction of a Cartesian coordinate system are assumed 
to allow only single-lane east-bound traffic while all those 
parallel to the ^-direction allow only single-lane north- 
bound traffic. Each of the lattice sites represents the 
crossing of a east- west street and a north-south street. 
In the initial state of the system, N x (N y ) vehicles are 
distributed among the east-bound (north-bound) streets. 
The states of east-bound vehicles are updated in par- 
allel at every odd discrete time step whereas those of 
the north-bound vehicles are updated in parallel at every 
even discrete time step following a rule which is a simple 
extension of the fully asymmetric simple exclusion pro- 
cess 0: a vehicle moves forward by one lattice spacing 
if and only if the site in front is empty, otherwise the 
vehicle does not move at that time step. Jamming arises 
from the mutual blocking of the flows of east-bound and 
north-bound traffic at various different crossings. The 
BML model has been modified and extended [flO| — |l7|| . 

We model the network of the streets as a N x N square 
lattice. The streets parallel to X and Y axes allow only 
east-bound and north-bound traffic, respectively, as in 
the original formulation of the BML model. A signal is 
installed at every site of this N x N square lattice where 
each of the sites represents a crossing of two mutually 
perpendicular streets. The separation between any two 
successive crossings on every street is assumed to consist 
of D cells so that the total number of cells on every street 
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is L = N x D. The linear size of each cell may be inter- 
preted as the typical length of a car; each of these cells 
can be either empty or occupied by at most one single 
vehicle at a time. Because of these cells, the network of 
the streets can be viewed as a decorated lattice. However, 
unlike the BML model || and the model of Horiguchi and 
Sakakibara |l4j which correspond to D = 1 and D = 2, 
respectively, D (< L) in our model is to be treated as a 
parameter. Note that D introduces a new length scale 
into the problem. 

The signals are synchronized in such a way that all 
the signals remain green for the east-bound vehicles (and 
simultaneously, red for the north-bound vehicles) for a 
time interval T and then, simultaneously, all the signals 
turn red for the east-bound vehicles (and, green for the 
north-bound vehicles). Thus, the parameter T introduces 
a new time scale into the problem. 

As in the original BML model, no turning of the vehi- 
cles is allowed. Therefore, the total number of vehicles 
on each street is determined by the initial condition, and 
does not change with time because of the periodic bound- 
ary conditions. 

Following the prescription of the NS model, we allow 
the speed V of each vehicle to take one of the V max + 
1 integer values V = 0, 1, V max . Suppose, V n is the 
speed of the n-th vehicle at time t while moving cither 
towards east or towards north. At each discrete time step 
t — ► t + 1, the arrangement of N vehicles is updated in 
parallel according to the following " rules" : 
Step 1: Acceleration. If, V n < V max , the speed of the 
n-th vehicle is increased by one, i.e., V n — > V n + 1. 
Step 2: Deceleration (due to other vehicles or signal). 
Suppose, d n is the gap in between the n-th vehicle and 
the vehicle in front of it, and s n is the distance between 
the same n-th vehicle and the closest signal in front of it. 
Case I: The signal is red for the car under consideration: 
If min(d n , s n ) < V n then V n — > min{d ni s n ) — 1. 
Case IT. The signal is green for the vehicle under con- 
sideration: 

There are two possibilities in this case: (i) When d n < s„, 
then V n — * d n — 1 if d n < V n , The motivation for this 
choice comes from the fact that, when d n < s n , the hin- 
drance effect comes from the leading vehicle, 
(ii) When, d n > s n , then V n — > min(V n ,d n — 1) if 
min(y n , d n — 1) x r > s„, where r is the number of the 
remaining time steps before the signal turns red. The mo- 
tivation for this choice comes from the fact that, when 
d n > s n , the speed of the n-th vehicle at the next time 
step depends on whether or not the vehicle can cross the 
crossing in front before the signal for it turns red. 
Step 3: Randomization. If V n > 0, the speed of the n- 
th vehicle is decreased randomly by unity (i.e., V n — * 
Vn — 1) with probability p (0 < p < 1); p, the random 
deceleration probability, is identical for all the vehicles 
and does not change during the updating. 
Step 4 : Vehicle movement. Each vehicle is moved for- 
ward so that for the east-bound vehicles, X n — » X n + V n 
where X n denotes the position of the n-th vehicle at time 



t while for the north-bound vehicles, Y n — > Y n + V n where 
Y n denotes the position of the n-th vehicle at time t. 

These rules are not merely a combination of the rules 
proposed by BML || and those introduced by NS || but 
also involves some modifications. For example, unlike 
all the earlier BML-type models, a vehicle approaching 
a crossing can keep moving, even when the signal is red, 
until it reaches a site immediately in front of which there 
is either a halting vehicle or a crossing. Moreover, if 
p = every east-bound (north-bound) vehicle can adjust 
speed in the deceleration stage so as not to block the 
north-bound (east-bound) traffic when the signal is red 
for the east-bound (north-bound) vehicles. 

In our computer simulations, we begin with an ini- 
tial configuration where N x and N y vehicles are put at 
random positions on the east-bound and north-bound 
streets, respectively. The states of the vehicles are up- 
dated in parallel following the rules mentioned above. 
After the initial transients die down, at every time step, 
we compute the average speeds (V x ) and (V y ) which are 
merely the averages of the instantaneous speeds of the 
east-bound and north-bound vehicles, respectively. The 
density c = (N x + N y )/(2LN - N 2 ) of the vehicles is the 
ratio of the total number of cars and the total number 
of cells in the system. Here we present the data for only 
the symmetric case N x = N y , for only a few set of values 
of the parameters £>,T, c,p, L, V max ; more details will be 
published elsewhere fig] . 

In the "free-flowing" phase of the BML model, both 
(V x ) and (V y ) oscillate between zero and a non-zero value 
periodically at odd and even time steps. In sharp con- 
trast, the time-dependences of (V x ) and (V y ) are much 
more realistic in our model, as is evident from fig. [I]. 
Moreover, as expected, V g , the maximum allowed val- 
ues of (V x ) and (V^) in the corresponding green phase, 
is smaller when the density c is higher. In this param- 
eter regime, following the switching of the red (green) 
signal to green (red), (V x ) rises (falls) to reach V g (0); 
the corresponding relaxation time is denoted by t g (t r ). 
For a given c, we now derive approximate analytical ex- 
pressions for t g and t r in terms of Vns{c), the steady 
speed of the vehicles, for the vehicle density c, in the NS 
model with periodic boundary conditions. Then, using 
the numerical estimates of Vns(c) from computer sim- 
ulations of the NS model we compute (V x ) and (V y ) for 
our model and compare with the numerical data obtained 
from direct computer simulation. 

We assume that during the red phase compact (i.e. 
without 'holes') queues of length N q — cD are formed in 
front of each signal. We now estimate the time t g until 
the stationary speed V g = Vns{c) is reached. There are 
two different additive contributions to t g . First, the last 
vehicle in a compact queue of N q vehicles starts moving 
after t\ = N q /(1 — p) timesteps since the leading vehicle 
in the remainder of the queue moves with probability 
I— p. Second, a halting vehicle reaches the speed Vivs(c) 
after a time t^ = Vjvs(c)/(1 — p) since it accelerates in 
each timestep with probability 1 — p. Thus, 
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FIG. 1. Time-dependence of average speeds of vehicles. 
The symbols o, +, □ and x correspond, respectively, to 
(14) for c = 0.1, (V y ) for c = 0.1, (14) for c = 0.5 
and (V y ) for c = 0.5. The common parameters are 
V max = 5,p = 0.1, D = 100 and T = 100. The continuous 
line has been obtained from the equations (^|) and (^). 
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where we have assumed that the green phase starts at t = 
0. Moreover, under the assumption that the relaxation 
is perfectly linear we obtain the following average speeds 
during a green phase: 



v g (t) 



f V NS (c)t/t g 
1 V NS (c) 



for t < tg 
for t > t„ 



(2) 



In the stationary state of the green phase, there are, 
on the average, - — I empty cells in front of each vehicle. 
So if the leading vehicle of a pair happens to be the last 
member of a queue already formed in front, then the 
following vehicle of that pair will move for to = (- — 
1)/Vns(c) with velocity Vjvs(c) and then stop suddenly 
since it reaches the tail of a queue. Therefore, on the 
average, it takes a time 



t r = cDto = 



V NS (c) ' 



(3) 



to form a queue of length N q — cD after the first vehicle 
has stopped at the red signal. In general, the vehicle 
nearest to a signal will not stop immediately after the 
signal turns red, but will keep moving for some time, 
say, t p . The distance of this vehicle from the signal is 
expected to be a fraction a of the average distance - — 1 
to the next vehicle ahead of it. Taking a = 1/2, for 
example, one obtains t p = 2cVns(c) ' Hence, the average 
speed during a red phase starting at t = is given by 



r v NS (c) 

Vr(t) = I V NS (C) 

1° 



t-tp 



for < t < t p , 

for t p <t <t p + t r , (4) 

for t > t p + t r . 



For the validity of these estimates T(2> 1) should be 
sufficiently large to guarantee complete queueing of the 
vehicles during the red phase, p should be small enough to 



ensure compactness of the queues, c should be sufficiently 
small so that the vehicles emerging from a queue should 
not be hindered by the halting vehicles of another queue 
in front. Moreover, the smaller is V max the stronger is 
the deviation from the linear relaxation assumed above. 
Furthermore, we have assumed that all queues have the 
same length. 

In fig. [y, T = 100 is sufficiently long so that (V x ) and 
(V y ) relax to Vns(c) during the green phase of the corre- 
sponding signals and to zero during the red phase of the 
signals. In contrast, for the same values of the parame- 
ters D and c, the average speeds (V x ) and (V y ) do not 
get sufficient time to relax to zero during the red phase of 
the corresponding signal provided T is sufficiently small, 
e.g., T= 20 (see fig. g). 
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FIG. 2. Time dependence of average speeds of vehicles. 
The symbols o, +, □ and x correspond, respectively, to 
(V x ) for T = 20, (V y ) for T = 20, (14) for T = 100 
and (V y ) for T = 100. The common parameters are 
Vmax = 5,p = 0.1, D = 100 and c = 0.5. The data for T = 20 
are shown only upto 40 time steps to avoid overcrowding of 
data points. The lines are to serve as guides to the eye. 

The most dramatic result of our investigation is that at 
a sufficiently large density c*(D, T), which depends on D 
and T, a phase transition from the " free-flowing" dynam- 
ical phase to a completely jammed phase can take place 
in our "unified" model. In the jammed phase, the flow 
of east-bound vehicles is blocked by the north-bound ve- 
hicles and vice-versa; this gridlock phenomenon as well 
as the typical configurations in the jammed phase (see 
fig. ||) of our "unified" model are similar to those in the 
BML model. In spite of these apparent similarities, as 
we shall explain now, the mechanism that triggers jam- 
ming in our "unified" model is different from that in the 
BML model. It is obvious from the updating rules that 
if p = 0, i.e., if no random braking takes place, com- 
plete jamming is impossible in this model at any density 
c < 1. Therefore, a vehicle which is located at the cross- 
ing of two mutually perpendicular streets and whose in- 
stanteneous speed is V = 1 at the end of the deceleration 
stage (i.e., having at least one empty site in front of it) 
would vacate the crossing unless its speed is reduced to 
V = because of random braking. However, if such a 
halt at a crossing due to random braking takes place at 
the last time step before the signal for it turns red, it 
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would not only continue to block the perpendicular flow 
of traffic through the same crossing during the next T 
time steps but would also give rise to a queue of jammed 
vehicles in the perpendicular street passing through the 
same crossing. Therefore, when p ^ and the density is 
sufficiently high, the dynamical phase of " freely-flowing" 
traffic becomes unstable against the spontaneous forma- 
tion of jams and the entire traffic system self-organizes 
so as to reach the completely jammed state. 
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FIG. 3. A typical jammed configuration of the vehicles. 
The east-bound and north-bound vehicles are represented by 
the symbols — > and |, respectively. 

Moreover, for given D, the shorter is the time interval 
T the smaller is the c* (see fig. |]). Besides, the density 
corresponding to the maximum flux also shifts to smaller 
densities with the decrease of T. Furthermore, the maxi- 
mum throughput is a non-monotonic function of T in the 
" free- flowing" phase; this result may be of practical use 
in traffic engineering for maximizing the throughput. 

In this communication we have developed a "unified" 
model where the jams are created by the same stochastic 
process as in the NS model, but the transition to com- 
plete jamming and the jammed configurations are very 
similar to those in the BML model. We have also es- 
tablished that our model describes the time-dependence 
of the average speeds of the vehicles in more realistic 
manner than any of the earlier CA models of BML- 
type. Results of our ongoing investigations of the ef- 
fects of (a) turning of vehicles from east-bound (north- 
bound) streets to north-bound (east-bound) streets and 
(b) green-wave signalling on the flow and jamming will 
be reported elsewhere Q . 

We thank Ludger Santen for useful discussions as well 
as for help in producing figure || and Dietrich Stauffer 
for valuable comments on the manuscript. This work is 
supported in part by the SFB34I Aachen- Julich-Koln. 
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FIG. 4. The fundamental diagram. The symbols O, +, □, x 
and A correspond, respectively, to T = 100, 50, 20, 10, 4. The 
common parameters are V ma x = 5,p = 0.5, L — 100, and 
D = 20. 
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